Pion structure function at small x from DIS data 
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Production of leading neutrons in DIS is usually considered as a tool to measure the pion structure 
function at small x accessible at HERA. The main obstacle is the lack of reliable evaluations of the 
absorption corrections, which significantly suppress the cross section. We performed a parameter 
free calculation within the dipole approach and found the absorption corrections to be nearly as 
strong, as for neutron production in pp collisions. We also included the significant contribution of 
the iso-vector Reggeons with natural (p, 02) and unnatural (01, p-n cut) parity with parameters 
constrained by phenomenology. With a certain modeling for the pion-to-proton ratio of the structure 
functions we reached good agreement with data from the ZEUS and HI experiments, successfully 
reproducing the observed dependences on the fractional neutron momentum z, the photon virtuality 
Q 2 , and the transverse momentum transfer qr- 

PACS numbers: 13.60.-r, 13.60.Rj, 13.60.Hb, 14.40.Be 



I. INTRODUCTION 

Neutron production in deep-inelastic scattering (DIS) 
on a proton can serve as a sensitive tool to study the 
properties of the meson cloud of nucleons, because only 
iso-vector quantum numbers in the crossed channel are 
allowed. If neutrons are produced at forward rapidities 
with small transverse momenta, the contribution of large 
impact parameters of j*p collisions dominates, so one 
can probe light mesons in the proton wave function, in 
particular pions. In terms of the dispersion relation this 
means that in this kinematic region one gets close to the 
pion pole. 

Thus, one can treat leading neutron production in DIS 
as a method to measure the structure function of the 
pion, F^(x^, Q 2 ), as is illustrated in Fig. [I] Summing up 
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FIG. 1: Graphical representation of the cross section of in- 
clusive neutron production in hadron-proton collisions, in the 
fragmentation region of the proton. 

all final states A" at a fixed invariant mass Mx one arrives 
at the total hadron-pion cross section at cm. energy Mx ■ 
This cross section is a slowly varying function of Mx, 
what leads to an approximate Feynman scaling. The 
rapidity gap covered by the pion exchange correspond to 
the energy, which is much smaller than the total cm. 
collision energy squared 
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where sq is the scale factor, usually fixed at 1 GeV 2 ; and 
z = Pn/Pp is the fraction of the proton light-cone mo- 
mentum carried by the neutron. If z is large, it is close 
to Feynman xf- 

The pion exchange brings in the cross section the factor 
(1 — z )~ 2oi7 '( t \ where a n (t) is the pion Regge trajectory. 
This factor is independent of the collision energy, if z is 
fixed, so the pion exchange contribution does not vanish 
with energy. The smaller is the 4-momentum transfer 
squared t, the closer one approaches the pion pole in 
the dispersion relation, and the more important is the 
pion contribution. However, the smallest values of t are 
reached in the forward direction at z — > 1. The lat- 
ter condition leads to the dominance of other Reggeons 
which have higher intercepts. Indeed, the correspond- 
ing Regge factor (1 — z ) _2qjr (*) for p and an Reggeons 
is about 1/(1 — z) times larger than the one for pion. 
Although in general these Reggeons are suppressed by 
an order of magnitude compared to the pion they 
become important at z > 0.9. 

The effective contribution of the axial- vector states (a\ 
pole and p-7r cut) with the parameters fixed from phe- 
nomenology, was found recently [2[ to be crucial explain- 
ing data on azimuthal asymmetry of leading neutrons 
produced in pp collisions. This Reggeon having a low 
intercept, affects the cross section at small z. 

The most important correction, which is the main fo- 
cus of this paper, is the effect of absorption, or ini- 
tial/final state interactions. The active projectile par- 
tons participating in the reaction, as well as the specta- 
tor ones, can interact inelastically with the proton target 
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or with the recoil neutron, and initiate particle produc- 
tion, which usually leads to a substantial reduction of the 
fractional neutron momentum. The probability that this 
does not happen, called sometimes survival probability 
of a large rapidity gap, leads to a suppression of leading 
neutrons produced at large z, because this process is asso- 
ciated with formation of a rapidity gap Ay ~ — ln(l — z). 
Some calculations predict quite a mild effect, of about 
10% even in the soft process pp — > nX while oth- 

ers |7Hld expect a strong reduction by about a factor of 
2. See [8] for a discussion of the current controversies in 
data and theory, for leading neutron production. 

Notice that usually the absorptive corrections are cal- 
culated in a probabilistic way, convolving the gap survival 
probability with the cross section. We found, however, 
that the spin amplitudes of neutron production acquire 
quite different suppression factors 0, lHj, and one should 
work with amplitudes, rather than with probabilities. 

At first glance the absorptive corrections of the 
hadronic fluctuations of a highly virtual photon should 
be vanishingly small. However, the observed weak Q 2 de- 
pendence of nuclear shadowing in DIS demonstrates that 
this is not true: both shadowing and absorption are dom- 
inated by rare soft fluctuations of the photon [ll|. This 
is why the absorptive corrections were calculated in Q 
relying on the effective absorption cross section adjusted 
to data on nuclear shadowing [12j |. 

Even more simplified evaluation of absorptive correc- 
tions were performed in p| H|, basing on the two com- 
ponent model for fluctuation of the virtual photon, soft 
and hard. The former was assumed to interact like the p 
meson, while the latter cross section was fixed zero. 

Below we perform explicit calculations of the absorp- 
tive corrections caused by the interactions of the fluctu- 
ations of a highly virtual photon within the dipole ap- 
proach citezkl. Moreover, like in pp collisions [10| . even 
a stronger absorption, related to the formation of a large 
color octet dipole in 7*7r interaction, affects the large-z 
part of the neutron spectrum. This effect has been missed 
in previous calculations of the absorption corrections. 

Below our results are presented as follows. In Sect. |TT] 
the spin structure of the amplitude without absorption 
corrections is presented, and the theoretical uncertain- 
ties in the evaluation of the cross section are discussed. 
Although the goal of the present paper is to study the 
possibilities of extraction of the pion structure function 
from data, we try to predict the cross section of leading 
neutron production, modeling the ratio of pion to proton 
structure functions. 

Sect. lIIIl is devoted to the absorptive corrections, which 
are the main focus of this paper. The important obser- 
vation is the production of a large size color octet dipole 
formed by the remnants of the virtual photon and pion. 
Initial/final state interaction of such a dipole, controlling 
the absorptive corrections at large z, only slightly depend 
on on the size of the qq fluctuation of the virtual photon, 
therefore almost no Q 2 is predicted. An examples of the 
absorption suppression factor Si q (b) and corrected for 



gluon radiation Sn q (b) are shown in Fig. |4j The same fig- 
ure demonstrates the reduced suppression factor 5 7 *(6) 
at smaller z, where absorption of the hadronic fluctua- 
tions of the virtual photon plays major role. A sizable Q 2 
dependence is predicted at smaller z, which is confirmed 
by data. The cross section of leading neutron production 
is found to be about twice smaller than the absorption 
uncorrected one, as is demonstrated in Fig. [7] 

The isovector Reggeons, which also contribute to the 
neutron production, are evaluated in Sect. IPV1 The high- 
intercept p-Reggeon is important at large z and large 
momentum transfer qx (it flips helicity). The low in- 
tercept ai-Reggeon contributes at smaller z. The Regge 
ai-pole itself is found to be very weak, and is replaced by 
an effective pole a\, which also represents the p-n Regge 
cut. 

All the contributions together, corrected for absorp- 
tion, reproduce data pretty well, as is demonstrated for 
the qr- and Q 2 -integrated cross sections in Fig. [5] The 
Q 2 and qx dependences are also well reproduced, as is 
shown in Fig. ® and [TU] respectively. 

The possibility of extraction of the pion structure func- 
tion from data is discussed in Sect. I VII and the sensitivity 
of the final results to the value of the pion-to-proton ratio 
of F 2 is demonstrated in Fig. [T3] 



II. BORN APPROXIMATION 

A. Measuring the pion structure function 

In the Born approximation the pion exchange contribu- 
tion to the amplitude of neutron production j*p —¥ Xn, 
depicted in Fig.[T] in the leading order in small parameter 
ttin/\/s has the form 



Ap-t n {q, z) — 



0-3 Ql 



-.a ■ q T 



£ p 4 B (qT,z), (2) 



where a are Pauli matrices; £ Pj „ are the proton or neutron 
spinors; qr is the transverse component of the momentum 
transfer; 



1-z 
Ql = — r=~ m N . 



(3) 



At small 1 — z <C 1 the pseudoscalar amplitude 
i> B (qT,z) has the triple- Regge form, 



x A r ^ x (M 2 x ), (4) 

where M\ = (1 — z)s, and the 4-momentum transfer 
squared t has the form, 



t = q\ + - Qt ■ 



(5) 
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and rjn(t) is the phase (signature) factor which can be 
expanded near the pion pole as, 



rfo(t) =i-ctg 



1 



(6) 



We assume a linear pion Regge trajectory a^t) = a' n (t — 
m 2 ) with a'^ ~ 0.9 GeV~ 2 . The imaginary part in ([6]) is 
neglected in what follows, because its contribution near 
the pion pole is small. 

The effective vertex function G 7r + p „(t) = 
g v+pn exp(R 2 t), where 4 +p „(*)/87r = 13.85. The 
value of the slope parameter R\ is specified below. 

The amplitudes in (J2])-(j4|) are normalized as, 



M 2 X ) 



(Mj<)| 2 , (7) 



where different hadronic final states X are summed at 
fixed invariant mass Mx- Correspondingly, the differen- 
tial cross section of inclusive neutron production reads 

mm. 



da 



p— fn 



dz dq 2 ^ 



(^) 2 |*|G^ pn (t)|r ?7r (t)| 2 (l-^) 1 - 2 -W 



x a t V (*' = Jl4). 



(8) 



The virtual photoabsorption cross section can be ex- 
pressed in terms of the structure function, 
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where 
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and x = Q 2 /s. 

Thus, the process of leading neutron production in DIS 
described by Eq. ^ , offers a unique opportunity to mea- 
sure the pion structure function at small x. 

Experimental data are usually presented in the form 
of ratio of neutron production and inclusive DIS cross 
sections flU Il7| , which in the Born approximation can 
be represented as 
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The last factor in the right-hand side of Eq. (fTTj) is 
known and provides a sizable suppression. Indeed, at 



small x < 0.01 the measured proton structure functions 
can be parametrized as 
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where A(Q 2 _) = a ln(Q 2 /A 2 ) with a = 0.048 and A = 
0.29GeV]l!]. F%( X7I ,Q 2 )/F%(x,Q 2 ) acquires a consid- 
erable suppression factor (1 — z) x ^ \ For example, at 
z = 0.8 and Q 2 = 13GeV 2 (the mean value in [l7[) this 
factor is 0.68. This factor is the main source of Q 2 de- 
pendence of the fractional cross section (fTTj) . which turns 
out to be pretty weak. For further calculations we rely 
on the more realistic QCD fit [19i ]. 



B. What to expect? 

The main unknown quantity in (jlip . which also is the 
main goal of experimental studies of this process, is the 
pion structure function, which enters the ratio (fl"2|) . Here 
we attempt at specifying the expected value of the ratio 

R^/ N (Xrr,Q 2 ), Eq. (US]). 

The hadron structure function F^x.Q 2 ) is pro- 
portional to the total cross section of interaction of 
the virtual photon with the hadron, F^ix^Q 2 ) = 

Q 2 1 (4Tr 2 a em ) <T] f(x, Q 2 ). In the target rest frame inter- 
action with a highly virtual photon is a perfect counter 
of the number of quarks in the target. Indeed, the in- 
teraction radius of a small (~ 1/Q 2 ) dipole is also small 
(~ l/lnQ 2 ), therefore interaction of the dipole simul- 
taneously with two target valence quarks, separated by 
a large distance, is suppressed. So the small dipole in- 
teracts separately with each target quark via a colorless 
exchange (Pomeron), i.e. the dipole-quark cross section 
is finite and universal, and the total dipole-hadron cross 
section is proportional to the number of the quarks. One 
arrives at the additive quark model, which was first pro- 
posed, though ill justified, for soft hadronic interactions 
20]. However, a highly virtual photon interacting with a 
large light hadron this model should be rather accurate, 
so one should expect R^/n = ^2/^2 = N^/N?, where 
Nfj is the number of quarks in the hadron h. 

One can also interpret this via the QCD evolution at 
small x. The cross section of interaction of a small-size 
qq dipole with a proton is proportional to gluon density 



crq q {r T ,x) = —a s {Q 2 )xg(x,Q 2 )rT, 



(14) 



where tt ~ 1/Q is the transverse dipole separation. 
There are many experimental evidences for existence in 
the proton a semi-hard scale, the mean gluon trans- 
verse momentum, of an effective gluon mass, of the or- 
der of Q ~ 700 MeV [HHIIj]. This means that glu- 
ons are located within a small distance ~ ?'o = 0.3 fm 
around the sources. Probing the proton at the semi-soft 
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scale Qq one resolves only the "constituent" quarks, but 
not their structure. At a higher scale the gluons and 
sea quarks are resolved as well, but the QCD evolution 
leaves them essentially within the same spots around the 
sources. Although the radius of the spots rises with 1/x 
as (r 2 ) = 4a'ln(l/x), the effective slope at a hard scale 
is small a' sa 0.1 GeV" 2 0-111], and the spots in the 
proton do not overlap up to the energy of LHC 0, HI| ■ 
Therefore, it is reasonable to expect that the amount of 
glue and sea quarks generated at small x though the evo- 
lution, is proportional to the number of the quarks, which 
are resolved at the soft scale, and play role of the initial 
condition for the evolution. 



1. 3 valence quarks in the proton 

In the non-relativistic quark model one may expect a 
simple relation, 



Eq. (jT5|) should be increased: 



Q) 



3' 



(15) 



provided that x is sufficiently small, and Q 2 is large. This 
value was used in all previous calculations of the cross 
section of leading neutron production in DIS. 

It is clear, however, that the relation (fT5)) is a simpli- 
fication, which misses the possibility of interaction with 
those constituents of the proton, which are different from 
just three valence quarks. Indeed, even the process under 
consideration is an example: as is depicted in Fig. [TJ the 
virtual photon probes quarks and antiquarks in the pion 
cloud of the proton. 



2. A multiquark proton 

A proton experiences quantum fluctuations to the 
states containing more than 3 quarks. This is pretty 
obvious at a hard scale, since the flavor-symmetric sea 
of quarks and antiquarks is generated perturbatively 
through the QCD evolution. Such a source of extra 
quarks ceases, at a soft scale, so one might think that 
the gluon density at small x in the proton at the starting 
semi-hard scale is proportional to the number of valence 
quarks, like is assumed in Eq. ([T5j) . 

There are, however, nonperturbative quantum fluctua- 
tions in the proton, which produce extra quarks at a soft 
scale, also contributing to the initial conditions for the 
evolution. One of such mechanisms is directly related to 
the process under consideration. Production of leading 
neutrons is a part of the inclusive DIS cross section and 
also contributes to F%. On the other hand, as one can 
see in Fig. [TJ the small dipole {qq} 7 * does not interact 
with the 3-quark nucleons via gluonic exchanges, but in- 
teracts with the pion, i.e. with a pair of extra quarks in 
the proton. Within the pion cloud model of the proton 
the number of quarks in the denominator of the ratio 



N p 
g 



3 + 2K), 



(16) 



where (n„) is the mean number of pions in the proton. 
The bottom bound for this correction is easy to estimate 
integrating the fractional cross section Eq. (fl2|) . 



0.15. 



(17) 



Since neutral pion exchange should also provide a half of 
this contribution, we can estimate, 



> 0.225. 



(18) 



This is the bottom bound because other final states, like 
baryon resonances, should also be added. This estimate 
is compatible with evaluations in [2_7j , which ranges from 
(tItt) = 0.25 to 0.38, dependent on the used model for the 
pion flux (uncorrected for absorption) and with earlier 
estimates in [28l - l30j . 

The mean number of pions can also be evaluated basing 
on the observed deviation from the Gottfried sum rule 
3l| of the measured flavor asymmetry of the proton sea, 



I as = j dx [d p (x) - u p {x)] (19) 



The E866 experiment at Fermilab measured the value of 
asymmetry at Ias = 0.118 ± 0.012 (32j, which results 
in the mean number of pions (n f ) = 0.36 [33], with the 
usual assumption that the weight of the |ttA) Fock state 
is half of that for the |7rA^) component. A somewhat 
larger values of flavor asymmetry, but with larger errors, 
were found in the NMC experiment, I A s = 0.148 ±0.039 
[H, and HERMES, I as = 0.16 ±0.03 experiments. 
The deduced expectations for the number of pions are 
(n OT ) = 0.44 and 0.48 respectively. 

The number of quarks in the proton gets contribution 
not only from the flavor asymmetric, like in Eq. (|16[) . but 
also from the flavor symmetric sea. At importance of the 
latter indicate data [32| on d v (x)/uy(x), which cannot be 
described by the pion cloud model |29l l33j , and need iso- 
scalar contributions, like a and ui mesons. The analysis 
of this data performed in (36| within the meson cloud 
model, conclude that data on d p {x)/u p {x) ratio can be 
described with the weight factors (n a ) = 0.023 — 0.078 
and (n u ) = 0.063 — 0.671, which may be rather large, but 
are quite uncertain. 

Inclusion of the flavor symmetric sea into the relation 
(|16[) might considerably increase the mean number of the 
proton constituents at the soft scale, 



ATP = 3 + 2«n 7r ) + <n (T ) + <7i w )). 



(20) 



Although the contribution of the isoscalar mesons may be 
significant, its magnitude is model dependent and poorly 
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known. Considering the above mentioned result of the 
E866 experiment, (n n ) = 0.36 as a lower value, we fix the 
total meson contribution at (n meson ) = 0.5 and and use 
it in the following calculations. With this value = 4 
and instead of the simplified expectation Eq. (TTSl) , we will 
rely on 



Rtt/n( x it,Q 2 ) — ~, 



(21) 



Although this number has a large uncertainty, we will 
rely on it through all further calculations up to compar- 
ison with DIS data for neutron production. Eventually 
we will check the sensitivity of data to Rk/n- 

Notice that we have not touched so far the nomina- 
tor of this ratio, the number of quarks in the pion, N„. 
Apparently, it might be also subject to corrections due 
to soft multiquark fluctuation, e.g. tt — > Trp, which is 
probably the strongest Fock component (it — > 2ir is for- 
bidden). However, pion is the Goldstone meson with an 
abnormally small mass, so any fluctuation is strongly 
suppressed by the energy denominator. In particular, 
the amplitude of the ir — ¥ np transition is suppressed as 
m^/(TOp + m w ) 2 . Therefore the weight factor for the \irp) 
Fock component is so small that can be safely neglected. 



3. More uncertaities 

Another theoretical uncertainty in Eq. (fTTj) is related 
to the slope parameter R\ of the pionic formfactor of 
the nucleon. It has to be fixed by phenomenology, but 
the results of model dependent analyses are quite di- 
verse [3rl - f4ll | and vary from zero to R\ = 2 GeV -2 . 
This uncertainty affects the magnitude of the fractional 
cross section Eq. (fTTj). especially at medium values of z. 
The forward cross section (l/<Ji nc )da B /dzdq^\ qT= o cal- 
culated in the Born approximation Eq. (|lll) is depicted 
in Fig. [2] by the strip between upper (R\ — 0) and bot- 
tom (R\ = 2GeV" 2 ) curves. The calculations are done 
at Q 2 = 14GeV 2 , which is the mean value for the DIS 
data [13 at Q 2 > 2 GeV 2 also depicted Fig. H As we 
mentioned above, the Q 2 dependence of the fractional 
cross section is quite weak. 

We see that even within the uncertainty in the pa- 
rameter Ri, the Born approximation significantly over- 
estimates data [Tj} at large z, where the pion pole is 
expected to dominate. For further calculations we fix 
R\ = 0.3 GeV" 2 , which was chosen in 0, i, i, H El| as 
most reliable. 



III. ABSORPTIVE CORRECTIONS 

Calculation of absorptive corrections, or initial/final 
state interactions, is quite complicated in momentum 
representation, where they require multi-loop integra- 
tions. However, these corrections factorize in impact pa- 
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FIG. 2: (Color online) The forward fractional cross section of 
neutron production in DIS calculated in the Born approxima- 
tion, Eq. |QTJ, with Rf = (upper curve) and Rf = 2 GeV" 2 
(bottom curve). Data points are the results of the ZEUS ex- 
periment [ItJ. 



rameters, 



fp-+n(b, z) = fp_+ n (b, z) Sib, z, s) , 



(22) 



where S(b, z, s) is the suppression factor caused by ab- 
sorption. Then one can Fourier transform the amplitude 
back to momentum representation, and the calculations 
are greatly simplified. So, we should first perform Fourier 
transformation of the amplitude Eq. ([2) to impact pa- 
rameter representation. 



A. Impact parameter representation 

The partial Born amplitude at impact parameter b, 
corresponding to has the form, 



CT 3 q L e^{b,z)-i^-ef{b, z ) 

Jzb 



where 



SJ>) 

(23) 



9 B (b,z) = / d 2 q T e^4 B ( qT ,z) 
N(z) 



l-/3 2 e 2 



[K (eb)-K (b/P)]; (24) 



?f(M) = l I d 2 q T e^(b-q)<p B (q T ,z) 



N(z) 
1 - ^ 2 e 2 



€ K 1 (eb)--K 1 (b/j3) 



(25) 



Here 



N(z) = \g^ pn z{\ - z)<^ + ^e- R ^A y .^ x (M%) 

(26) 



G 



Z (<il + m l) ! 

- [Rl-a'v Hl-z)] 



(27) 



To simplify the calculations we replaced here the Gaus- 
sian form factor, exp(— f3 2 q 2 ^), by the monopole form 
1/(1 + f3 2 q^), which is a good approximation at the small 
values of qx we are interested in (both shapes are ad hoc 
anyway). At the same time we retain the Gaussian de- 
pendence on gi, which can be rather large. 



B. Survival amplitude of a {qq}1 -{qq}& dipole 

At large z — > 1 the process under consideration is 
associated with the creation of a rapidity gap, Ay = 
| ln(l — z)\, in which no particles are produced. Absorp- 
tive corrections, caused by initial and final state interac- 
tions of the projectile partons with the target and recoil 
neutron, may substantially reduce the probability of gap 
formation. Indeed, any inelastic interaction (color ex- 
change) of the active or spectator partons should cause 
intensive multiparticle production filling the gap. Usu- 
ally the corrected cross section is calculated probabilis- 
tically, i.e. convoluting the cross section with the sur- 
vival probability factor (see Q and references therein). 
This recipe may work sometimes as an approximation, 
but only for ^-integrated cross section. Otherwise one 
should rely on a survival amplitude, rather than probabil- 
ity. Besides, the absorptive corrections should be calcu- 
lated differently for the spin-flip and non-flip amplitudes 
(see below). 

The DIS on a virtual pion shown in Fig. [TJ i.e. the 
inelastic collision 7* + n —> X, can be seen as a color 
exchange between the colorless qq Fock component of 
the proton and the pion mediated by gluonic exchanges. 
Nonperturbatively, e.g. in the string model, the hadron 
collision looks like intersection and flip of strings. The 
final state of such a collision is two color octet qq pairs, 
originated from the photon and pion respectively, as 
is depicted in Fig. [3] Hadronization of the color-octet 
dipole, {qq}s — {qq}s, leads to the production of different 
final states X. 

According to Fig. [5b the produced color octet-octet 
state can experience final state interactions with the re- 
coil neutron. On the other hand, at high energies multi- 
ple interactions become coherent, and one cannot specify 
at which point the color-exchange interaction happens, 
i.e. initial and final state interactions cannot be disen- 
tangled. In terms of the Fock state decomposition the 
projectile proton fluctuates into a 4-quark color octet- 
octet before the interaction with the target. The fluctua- 
tion life-time, or coherence time, rises with energy and at 
high energies considerably exceeds the longitudinal size 
of target proton (see, however, more detailed discussion 
below). 

This leads to a different space-time picture of the pro- 



cess at high energies, namely: long in advance the in- 
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FIG. 3: (Color online) a: Born graph with single pion ex- 
change and excitation of the projectile photon, 7* + n — > X; 
b: inelastic interaction, 7* + n — > X, via color exchange, 
leading to the production of two color-octet qq dipoles which 
hadronize further to X; c: Fock state representation of the 
previous mechanism. A color octet-octet dipole as a 4-quark 
Fock component of the projectile photon, 7* — > {qq}s — {qq}$, 
interacts with the target proton via ir + exchange. This 4- 
quark state may experience initial and final state interaction 
via vacuum quantum number (Pomeron) exchange with the 
nucleons (ladder- like strips). 



teraction the incoming photon fluctuates into a 4-quark 
state 7* — > {qq}8-{<jq}8 which interacts with the target 
via pion exchange, as is illustrated in Fig. (3J:. The sur- 
vival probability amplitude Si q (b, r, s) for a dipole of sep- 
aration r colliding with a nucleon at impact parameter b 
can be estimated on analogy with [l(| as, 

S iq {b, r) = 1 - Im f 4q (b, r) « [l - Im f qq (b, r)] \ (28) 

Here we rely on the large N c approximation and re- 
placed a color octet-octet dipole by two triplet-antitriplet 
dipoles. 

Notice that the mean {qq}J -{qq}s separation is 
large, and interaction of the 4g-system with the nu- 
cleon target is soft, although the size of the {qq}J pair 
maybe as small as 1/Q- Therefore, for the dipole am- 
plitude f qq (b,f,s) we employ the parametrization de- 
pendent on energy, rather than Bjorken x. The b- 
integrated phenomenological dipole cross section Eq. (|52")l 
is parametrized in the saturated form [22| . 



a qq (r, S )=a (s) fl-e^W 



where R Q (s) = 0.88 fm (s /s) 014 ; s = 1000 GeV 2 ; 



o-o(s) 



(*) 1 + 



(29) 



(30) 



a tot( s ) = 23.6mb x (s/s ) 08 , and the mean pion charge 
radius squared is {r 2 h )^ = 0.44 fm 2 

The partial amplitude f qq (b, f, s, /3) of elastic scattering 
of a qq dipole with transverse separation r and fractional 
light-cone momenta (3 (for q) and 1 — /3 (for q) on a proton 
at impact parameter b was derived in [l^, EH, El j 
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Im fq q (b,r,s,P) = 



<7 (S 


) f 


(b-rP) 2 ~ 






SnB(s) \ 6XP 


2B(s) 


+ exp 






r 2 


[6+ (1/2 - 






2exp 


R 2 (s) 


2B(s) 


i 



[b + r(l-/3)f 
2B(s) 



(31) 



This partial amplitude satisfies the condition, 

2 / d 2 b Im fgJb, r, s, (3) = a 9 J r , s), (32) 



The fractions f3 and 1 — /3 of the light-cone momentum 
of the 4<7-system carries by the color octets {qq}J and 
{qq}g, which are the debris of of the photon and pion 
respectively, are related to z = 1 — M x j ' s as, 



Q 2 



8(1 



(33) 



where mj- is the transverse mass of the pion debris {qq}$ , 
which we fix at (m^) = 1 GeV 2 . 

The amplitude Eq. (|3"TT) also correctly reproduces the 
elastic ir-p slope B^(s) provided that the effective slope 
parameter B(s) has the form [loj . 



B(s) 



We use the Regge parametrization B^f(s) — B + 
2a'jp\n{s/fi 2 ), with B = 6GcV~ 2 , a' p = 0.25 GeV -2 , 
and p? = lGeV 2 . 

To get the differential cross section the absorp- 
tion corrected Born amplitude of neutron produc- 
tion, f p ^ n (b,f,z) = fp_+ n {b, z) S± q (b,r, z, s), should be 
Fourier transferred back to the momentum representa- 
tion, squared and averaged over the dipole size r, 



\fp-m(qT,z)\* 



d 2 rW 2 (r,M x ) 



dzdq 2 (2^ 
d 2 b d 2 b' e** & - P) f^ n (b, r, z)f p ^ n (b\ r,z), (35) 



where W 2 (r, M x ) is the probability distribution of im- 
pact parameter of 7*7r collision at cm. energy Mx- To 
simplify numerical calculation we employ here the same 
approximation as in [lol ] assuming that each amplitude 
can be averaged over r separately, i.e. the absorption 
factor Si q (b, z,s) in (|22|) is related to one in (f28| as 



(S 4g (&,w)) = / d 2 rW(r,Q 2 ,M 2 x )S iq (b,f, Z , S ). 



(36) 



To proceed further we have to specify the distribu- 
tion W(r, M x ) over the size r of the {qq}J -{qq}$ dipole, 
which is the impact parameter of the j*ir collision at cm. 
energy Mx- Therefore, the r-distribution W(r,M x ) 
is given by the partial elastic photon-pion amplitude 
P el v (r,M x ), for which we use the normalized Gaussian 
r-dependence, 



W(r,Q 2 ,M 2 x ) 



exp 



-r 2 /AB^{Q 2 ,M 2 



X ' 



(Q 



2 M 2 X ) 



(37) 



With a good precision, checked numerically 

(S 4q (b,r,z,s)) r w (l- lmfqq(b,r,z,s)) r , where 



(ka.fqq(b,r,z,s)) r 



Ms) J 1 

8tt I Bp(a,z) 



exp 



2Bp(s,z) 



1 



■ exp 



B 1 - (s,z) [ 2Bi~p{s,z)_ 
b 2 



B 6 {a,z) [l + ABy(M 2 x )/R 2 (s) 



■ exp 



2£ 5 (s,z) 



(38) 



where 



Bp(s,z) 



B(s) + 2f3 2 B2r(M 2 x ). 



(39) 



B 6 (z,s) = 5( S ) + 2£ 2 B^(M 2 ); 



l+ABy{M x )/B4{s) 



(40) 
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The Q 2 and energy dependences of the elastic 7*7r 
slope, B2 l 7 '(Mx,Q 2 ) are expected to be similar to 
what has been observed for electroproduction of dif- 
ferent vector mesons, p, <p, J/^, in j*p interactions 
[45M47I ] . It was found that the value of the slope sat- 



urates at Q 



"«v ~ 5GeV 2 at the universal level, 
2a' P ln(W 2 /fi 2 ), with B$ p a 4GeV" 2 
and a'p ~ 0.1 GcV~ 2 . Notice that the observed small 
value of the Pomeron trajectory slope is in a good ac- 
cord with the theoretically expectation and is consid- 
erably smaller than a' P rs 0.25 GeV~ 2 observed in soft 
processes, which are affected by saturation of unitarity 
[241 [25l | . For the 7*7r slope at high Q 2 we use a similar 
form 



Q 



2a' P \n(M 2 x / p 2 



(41) 



with Bl 71 



B 1 p 



2 GeV~ . The latter relation is is 
written in analogy to the systematics of slopes observed 
in soft hadronic collisions. 

The result of numerical calculation of S± q (b, z) at 
Vs = 100 GeV, z = 0.7 and Q 2 = 13GeV 2 are shown 
as function of b in Fig. U by the upper solid curve. The 
suppression factor hardly varies with z and is slightly 
enhanced with energy. Naturally, the absorption effect 




FIG. 4: (Color online) Partial survival amplitudes S4 q (b) (up- 
per solid curve) calculated with Eq. (|36l) and S l 4< ? (6) (bottom 
curve) calculated with Eq, (|44l) . at cm. energy of 7*-proton 
collision ^/s = 100 GeV and z = 0.7. Dashed curves show the 
suppression factor S^* (b) calculated for virtual photons with 
Q 2 = 3 (bottom) and 40 GeV 2 (upper dashed curve). 

is strongest for central collisions (suppression down to 
40%), and gradually ceases towards the periphery. 



Corrections for eluon radiation 



graph shown in Fig. [5] (left) this Fock state would corre- 
spond to the Pomeron gluonic latter without rungs, i.e. 
without gluon radiation. Gluon bremsstrahlung, is how- 
ever an important process, which is responsible for the 
observed rise of the cross sections with energy. So higher 
Fock states, like one depicted in Fig. [5] (right), contain- 
ing gluons should be taken into account. Apparently they 
should lead to enhanced absorption effects. 

In large- N c approximation a Fock component contain- 
ing gluons can be replaced by a multi-dipole state [48j . 
The Pomeron is known to have a two-sheet topology 
(cylinder) , which corresponds to the replacement of the 
4-quark state by two qq dipoles as was done above. Ev- 
ery gluon radiated within one of the two Pomeron sheets 
creates an extra color triplet 33 dipole. Correspondingly, 
the absorption factor gains an extra suppression factor 
1 — Im/g 9 (ro). Thus, a Fock state containing besides 
4-quark n gluons provide a modified absorption factor 

q(n) 

S iq (b) 



1 



(fqq{b,r gg )) 



(42) 



Such states containing gluons, are characterized by two 
different scales, two typical dimensions [23|, [24j. One is 
the large size of the pion, ~ 1/Aqcd, which dictates 
the mean size of the 4-quark color octet-octet dipole 
in (|37p . Another, much smaller distance is the mean 
size of a glue-glue dipole. Analysis of data on large- 
mass diffraction [22[ , as well as many other experimental 
observables [24[ show that this distance is quite small, 
(r gg ) = ro = 0.3 fm. Thus, averaging of the dipole ampli- 
tude, which has the same size r gg as the glue-glue dipole, 
should be averaged as, 

1 



d 2 r n 



^..rt f ?q (b,r gg ) (43) 



Further, we assume that the number of radiated gluons 
has the Poisson distribution, so we can sum up the ab- 
sorption factors of the s=Fock states with different num- 
ber of gluons, 



S iq (b) = S iq (b)e 



E 

n„=0 



n g \ 



1 - (fqq0,f gg )) 



= S iq (b)e-^U^ . 



(44) 



Jt, R.' 





So far our consideration was restricted to the low- 
est 4-quark Fock state of the photon, {gg} 7 » — {qq}-m 
contributing to neutron production. In the triple-Regge 



FIG. 5: (Color online) Left: triple-Regge graph nirlP for lead- 
ing neutron production in DIS; Right: the unitarity cut of the 
of this graph, shown as a gluon comb. 
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The mean number of radiated gluons (n g ) can be esti- 
mated looking at the X- dependence of the DIS cross sec- 
tion. In the leading-log approximation integration over 
rapidity of each gluon results in a factor ln(l/a;). Sum- 
ming over gluon number one gets for the total 7*-proton 
cross section, 



.7 v 



E 

n„=0 



[ffln(l/x)p 



(45) 



where g includes the coupling and other factors acquired 
due to radiation of each gluon. 

Thus, according to Eq. f| 13[) the mean number of radi- 
ated gluon reads, 



(n g ) = A(Q 2 )ln(l/^), 



(46) 



where A(Q 2 ) is defined in (|TT 

Now we are in a position to calculate the modifies ab- 
sorption factor Eq. (|44l) . The result is plotted by the 
bottom solid curve in Fig. [4] hi comparison with the un- 
corrected survival amplitude Si q . Although the modified 
suppression is stronger, as was anticipated, the difference 
is rathe small. This is a result of smallness of ro . 



D. Coherence lengths for the photon Fock states 

We have assumed so far that all the Fock components 
of the photon considered above have the lifetime, or co- 
herence length, much longer than the dimension of the 
target. This is certainly true for the simplest Fock state 



7* — > qq, which has a long coherence time, called Ioffe 
time (see more accurate evaluation in [Hj]), 



n q = 



1 



2xmpj ' 



(47) 



where x ~ 10 3 in the kinematics of HERA. 

However for the 4-quark Fock states this is not obvious. 
The coherence length is given by 



1 

<1L 



(1 - z)m N ' 



(48) 



This coherence time becomes very long for large z — > 1 

Since the target nucleon size is ~ 1 fm, only at 
z > 0.8 the coherence length Eq. (|4"5)l is sufficiently long 
to rely on the above evaluations of the absorptive effects. 

In another limit of a very short coherence length if 1 <C 
fjv not only initial, but also final state interactions of the 
4-quark state are impossible, because even the the inter- 
action time and time scale of creation of this state exceed 
the nucleon size. In this case the absorptive corrections 
are generated only by the interaction of long-living qq 
fluctuations of the photon. In this case the absorption 
factor has the form, 



S* 7 »(6) = 1 - (lmf qq (b,r,a)) 7 



(49) 



where the partial elastic amplitude Im f qq (b, f, a) is given 
by Eq. (|3"Tj). The averaging over the transverse dipole 
separation r, and the fractional light-cone momentum of 
the quark, a, is done as follows, 



(Im. f qq (b,r, a)) ra = 



da / d 2 ra qq (r,x 7r ) \^ qq {r,a,Q 2 ) 



da j d ra qq (r 1 x w ) \$! qq (r,a, Q 2 )\ luif qq (b,f,a), 

(50) 



where the weight factor l^ggl 2 = j^^l + l^ad. contains 
the standard photon distribution functions |50l . l5lj . No- 
tice the importance of the factor <j qq (r, x^ ) , which comes 
from the Born amplitude of neutron production. With- 
out this factor the result of (|50|) would be zero (49j . 
because the normalization of the distribution function 
of transversely polarized photons, 4^ is ultraviolet di- 
vergent. This divergency corresponds to ultra- heavy qq 
fluctuations, which dominate in a transversely polarized 
photon in vacuum. However, such fluctuation are "ster- 
ile" , i.e. cannot interact (color transparency) , while the 
process under consideration contains at least one dipole 
interaction. 

Examples of the results for S' 7 » (b) calculated at Q 2 = 3 
and 40GeV 2 are plotted in Fig. [4] by dashed curves. Ap- 
parently, the Q 2 dependence of the absorption factor is 



rather weak. This was anticipated, because in (|50p we av- 
erage the dipole amplitude squared, similar to diffraction 
of nuclear shadowing. The integral over a turns out to 
be dominated by the aligned-jet configurations [Hj], i.e. 
by the endpoint behavior of the distribution functions, 
a < m 2 /Q 2 [53l | . The corresponding qq transverse sep- 
aration becomes rather large, ~ l/m q , and independent 
of Q 2 , (r 2 ) - Q 2 a(l - a) + m 2 . We fixed the effec- 
tive quark mass, which is in fact the infrared cutoff, at 
m a = 0.15 GeV adjusted to data on nuclear shadowing 

Thus, we know the absorption corrections in two lim- 
iting regimes: (i) l\ q ^> rjv, in this case the suppression 
factor S4 q (b) is given by Eq. (|44| and depicted by solid 
curves m Fig.lU (ii) < r N , in this case the absorption 
factor S*y* (b) is given by Eq. (|4^| and is shown by solid 
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curves in Fig. |4j In order to interpolate between these 
limiting regimes we employ the following simple proce- 
dure, 

S(b) = S 4q (b) F N (q L ) + S r (b) [1 - F N (q L )] , (51) 

where the transition formfactor is chosen in the dipole 

-l 



The 



form, F N (q L ) = (l + q\L 2 ) 

parameter L characterizes the dimension of the target 
nucleon, so it should be of the order of 1 fm, and we fix it 
at this value, L = 1 fm, for further calculations. However, 
this parameter can be varied within a reasonable range. 
The z-dependence of the suppression factor Eq. (fSTj) is 
plotted in Fig. [5] for few values of b. 




FIG. 6: (Color online) The absorption factor S(b, z), Eq. J5TJ 
as function of z at fixed impact parameters b = 0, 0.5 and 
lfm. 



Notice that although the recipe Eq. (fSlj) is partially 
ad hoc, it interpolates between the two known limiting 
regimes of very long, l^ q 3> rpj {z > 0.8), and very short, 
<C rpj (z ~ 0.5) coherence length. As far as these two 
regimes are predicted, the interpolation procedure should 
not affect the results significantly. 



E. Cross section corrected for absorption 



As soon as the absorption factor Eq. (151]) is known, we 
can perform the inverse Fourier transformation to mo- 
mentum representation, 



fp->n(QT>z) = 



1 



d 2 be^ Tb 



(2tt)* 

x f^ n (b,z)S(b,z,Q\s) 



(52) 



where the Born amplitude in impact parameters, 
fp-> n (b,z) is given by ([23]) . So the absorption corrected 
partial spin amplitudes read, 



9 , a {b,z) = e§.(b,z)S{b,z). 



(53) 



Then the Fourier transformed amplitude reads, 

<Ps(qr,z) 



A p ^ n (q T ,z) = 



03<7l </>o(<?t, z ) + vqr- 



where according to (|2gj). pjj) and d^HJ) 

OO 

N(z) 



(54) 



27r(l-/3 2 e 2 ) 



K (eb)-K 



db b Jo(bqT) S(b, z) 



(55) 



qT& s {qT,z) 



N(z) 



dbb Ji(bq T )S(b,z) 



(56) 



Now we can calculate the differential cross section of 
inclusive production of neutrons corrected for absorption, 



dOr, 



dz dq\ 



ao(z, qr) + cr s (z, qr) , 



where 



oo{z,q T ) 



^ I0o(^,^)| 2 
s 

<T s (z,q T ) = — \cj) s (q T ,z)\ 2 
zs 



(57) 

(58) 
(59) 



The effects of absorptive corrections are illustrated in 
Fig. [7] The effect of suppression caused by the absorption 




0.6 0.7 0.8 0.9 



FIG. 7: (Color online) Forward fractional cross section of neu- 
tron production in DIS calculated without absorption correc- 
tions (dashed curve); absorption suppressed by factor Siqib) 
Eq. (|44|l (thin solid curve); and with the ^-dependent absorp- 
tion factor S(b,z) Eq. (|51[) (thick solid curve). 

factor S4 q (b) defined in (1581) is demonstrated by compar- 
ison of the fractional cross section of forward neutron 
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production (dashed curve) with absorption suppressed 
result plotted by the thin solid curve. We observe a 
rather strong effect, absorption reduces the cross section 
by nearly factor 2. Inclusion of the coherence time effect 
results in a z-dependent absorption factor S(b, z) defined 
in (j5Tj) and illustrated in Fig. [6j This final absorption 
corrected cross section is plotted in Fig. [7] by thick solid 
curve. We postpone comparison with data, because sev- 
eral more mechanisms of neutron production are to be 
added. 



IV. OTHER REGGEONS 

Besides pion exchange, other iso-vector Reggeons con- 
tribute to neutron production. Those are subdivided to 
natural parity Reggeons (p, 0.2)1 which have high inter- 
cepts Oe(0) rs 1/2, and unnatural parity ones (ai, n', 
etc.) with lower intercepts 



A. Natural parity Reggeons 

The leading Reggeons contributing to neutron produc- 
tion are p and 0-2. At large invariant masses M x they 
do not interfere with the pion exchange and with each 
other. Indeed, summing over final states at fixed Mx 
one gets the imaginary parts of the amplitudes for the 
processes, 7* + it —} 7* + p, or 7* + p — >• 7* + 0-2, which 
are suppressed by a power of 1/Mx- 

These amplitude are known to be dominated by their 
spin- flip part [54l [Hj], so we neglect the non-flip term 
in what follows. In the Born approximation the leading 
Reggeons contribute to the cross section as, 

1 (d^ n \ 2 /^' (l-^-M*) 
a mc \dzdq% ) r qT \ 8 J z 2 

x G^pn(«)M01 a ^ff • (60) 

We consider two leading exchange degenerate Reggeons 
p and ci2- The signature factor of the latter n a2 (t) — 
—i — cot [§Q! a3 (t)] diverges at the so called nonsense 
wrong signature point, to = — a a2 (0)/a' a ~ — 0.6GeV 2 . 
In order to kill this unphysical pole the residue function 
of the a,2 Reggeon must have a zero at this point, i.e. a 
factor (1 — t/to). According to exchange degeneracy the 
residue function of the p-Reggeon should also contain this 
factor, which is not compensated by any pole at t = to. 
This is confirmed by data on differential cross section of 
reaction n~p — > 7r°n, which indeed has a dip at t = to 

In this circumstances the i-dependences of p and 02 
exchange amplitudes are rather uncertain. Since we are 
focused on small-i region, the most reasonable solution 
seems to be to fix the signature factors of both Reggeons 



at rjn(0). Moreover, basing on the exchange degener- 
acy of p and 0,2 Reggeons, we fix a p (t) = a a2 (t) and 

G p+pn (t) = G a + pn {t). 

The contribution to the cross section of the spin-flip 
Reggeon amplitude can be described by the term similar 
to Eq. (l59|) . properly modified, 

-&«<«--'>i'- < 61 » 

where 4>%(qT,z), compared to Eq. (|56p , contains the 
imaginary part neglected for pions, and several modi- 
fications, 

00 

n (z) r 

<%(Qt,z) = / dbb JiMK^b/PJSfaz). 

2nq T pi J 


(62) 

The further notations are, 

N p (z) = ^g p+pnTlp (Q)z(l-z)- a ^ +a >l 

x e- R2 p^A Tp ^ x {M 2 x ) (63) 
0j, = \ K - ol p ln(l - z)\ . (64) 

The absorptive corrections for the Reggeons are calcu- 
lated with the same suppression factor S(b) in (|6"2"1) . Its 
effect maybe somewhat stronger for Reggeons than for 
pions, since the interaction in this case is more central. 

For the vertex function G P NN(t) = g P NN exp(R^t) we 
rely on the phenomenological global Regge analysis [55[ 
of high-energy hadronic data, which results in g P NN — 
0.5 g7rNN , and R 2 p = 1 GeV~ 2 . 

Notice that N p (z) contains an additional z- 
dependence, a factor ~ 1/yT — z, compared to the 
pion exchange, Eq. (|2T)1) . This mean that at z — > 1 the 
p/a,2 Reggeon contribution should dominate over the 
pion exchange [Jj . although at which z this happens 
depends on the couplings strength. 

B. ai-like exchanges 

The study of spin effect in leading neutron produc- 
tion performed recently |2J revealed important role of 
the axial-vector Reggeons, like a\ meson. Interference 
of the related non-flip spin amplitude with spin-flip pion 
exchange well explained data on transverse single-spin 
asymmetry of leading neutrons produced in pp collisions. 

If fact, the situation with axial vector mesons is more 
complicated. Assuming vector meson dominance in the 
axial current, in analogy with the vector current, one ar- 
rives at a dramatic contradiction of the Adler relation for 
diffractive neutrino-production of pins with data, the ef- 
fect called Piketti-Stodolsky puzzle (56[. It was proposed 
in [57| that the source of the problem is the assumed 
axial-vector dominance, while in reality the ai pole is a 



12 



very weak singularity, and the main contribution to the 
dispersion relation for the axial current comes from the 
p-ir cut. Indeed a detailed analysis of data on diffrac- 
tive dissociation n — > pn performed in Q shows that the 
invariant mass distribution of the produced in the l + 5 
wave p-n forms a pronounced narrow peak at a mass 
AIa 1 = 1.12 GeV close to the a\ mass. In many instances 
one can tread such a p-ir Regge cut as an effective a \ -pole 
0, [H^-IHH . Its contribution to the spin nin-flip part of the 
Born amplitude Eq. (0) reads, 



B _ a di 



Ga ipn {t) Vhl (t)(l - Z )- a ^ 



x A 7 . ai ^x(M|), 



where 



Vatit) =-i- tg 



Tra~ ai (t) 



(65) 



(66) 



The Regge trajectory of the p-ir cut has the form, 
a~ ai (t) = a„ P {t) = a n (0) + a p (0) - 1 + , * p , t, (67) 



so a 3l (Q) = -0.5; a' di = 0.45 GeV" 2 . 

The aiNN vertex is parametrized as Ga ipn (t) = 
oa lP „exp(i?| i t). The a^NN coupling was evaluated in 
Q basing on PCAC and the second Weinberg sum rule, 
in which the spectral functions of the vector and axial 
currents are represented by the p and the effective a\ 
poles respectively. This leads to the following relations 
between the couplings, 



2 p 

gtnNN _ m a 1 Jtt 
9ttNN 2mjv fp 



0.5, 



(68) 



where /„. = 0.93771^ is the pion decay coupling; f p — 
\f2m 2 p /^ p , and j p is the universal coupling (pNN, pnn, 
etc), 7 2 /4tt = 2.4. 

Applying to the Born amplitude Eq. (f65j) the proce- 
dure of correcting for absorption, developed above, we 
Fourier transform the Born amplitude to impact param- 
eters, introduce the absorption factor 5(6), the transform 
the result back to momentum representation and get, 



4>o{<lT,z) 



where 



dbbJ {bq T )K {b/p- ai )S{b,z), 



(69) 



B 2 



7r a 



x e- R ^ ql A Y ~ ai ^ x {M 2 x ) 
= ~ [ R l - 4x ln(l - z)] 



(70) 
(71) 



For the sake of simplicity we assume that A^&^x (Mjc) 
is equal to the same amplitude on pion, p, or a,2 targets. 



Besides, the interference between spin non-flip ampli- 
tudes with pion and a exchanges also contributes. Alto- 
gether the corresponding part of the cross section reads, 



*i)(*,«r) = f [\<t>l{q T ,z)\ 2 +e\4 1 {qT,z)\ 2 



2£Re<^ Ql (<z T ,2) 



(72) 



where £ = 2m x / \f\t\ ~ 2/(1 — z) is a factor related to 
the spin structure of the axial- vector vertex Q ■ 

In the interference term one needs to know the off- 
diagonal diffractive amplitude, 

X^W-fxA'*-** = M% ImA y . w _ Hy .g 1 (Af|- ) pr = 0) 
x 

(73) 

The amplitude of the process 7*7r — > j*di can be related 
to the reaction irp — > hip relying on Regge factorization, 
and dominance of the diffractive excitation 7r — > pir (see 
above). The amplitude Eq. ([75)) is suppressed compared 
with 7*7r elastic scattering by the factor f2 defined as, 



2 _ 



G?cr(7*7r — > j*ai)/dpj, 



da{"f*TT — > 7*7r)/dpj, 
da(np — » irpp)/dpj 



Pt=0 



dcrijrp Ttp)/dp7p 



Pt=0 



(74) 



According to L is energy independent and can be 
evaluated at M x — 150 GeV 2 where data for diffractive 
7rp interactions are available, do(jtp — > a\p) I^P%\ Vt —q — 

1.67mb/GeV 2 . Then, 



= 



16tt x 1.67 mb/ GeV 2 

1 Wp : 



= 0.29. 



(75) 



Thus, the third term in (f72"j) can be presented as, 

tan [7ra 5l (-g|)/2] 



(q T ,z) = -20 



\V~ai 



x yJ\^{q T ,z)Y\^{q T ,z)\ ,(76) 

where (t) is given by Eq. (|67]) and we neglected 
the small imaginary part of the pion signature function 
Eq. ©. 

Notice that the contribution of the a\ exchange to neu- 
tron production has been well tested. Analogous calcula- 
tions [IJ of the imaginary part of the it — &i interference 
led to a very good agreement with data on azimuthal 
asymmetry of neutrons produced by polarized protons. 



V. NUMERICAL RESULTS 

A. z-dependence 

Two experiments at HERA, ZEUS 0, [T^ and HI 
[60| have studied leading neutron production in DIS. 
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Their results are available for comparison with the same 
kinematics, integrated over within a fixed polar an- 
gle, <7t < z x 0.69 GeV, and up to a fixed maximum, 
qx < 0.2 GeV. These data are depicted by round points 
in Fig. [8] in the upper and bottom panels respectively. 
Analogous results of the HI experiment are published in 
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FIG. 8: Fractional differential cross section of neutron produc- 
tion, integrated up to q™ ax = 0.2 GeV (upper panel) and up 
to q^ ax = z x 0.69 GeV (bottom panel). ZEUS data Fl7| are 
shown by round points. Data from the HI experiment [601 ] de- 
picted by squares are normalized as described in the text. In 
all sets of data the error are dominated by the systematic un- 
certainties. The pion pole contribution, calculated including 
absorptioncorrections is plotted by long-dashed curve (com- 
pare with Fig.0). The contribution of p and a 2 Reggeons and 
the effective a\ pole are shown by short-dashed and dotted 
curves respectively. The sum of all these terms is presented 
by solid curve. 



the one measured by ZEUS at z < 0.8. This indicates at 
a significantly different gr-slopes of the differential cross 
section measured in these two experiments. 

The absorption corrected pion Regge pole contribu- 
tion, which was shown by thick solid curve in Fig. 
is depicted here by long-dashed curves. As was antici- 
pated, the contribution of the effective ai-Reggeon plot- 
ted by dotted curve, is vanishing at large z because of 
the low Regge intercept. On the contrary, the p and a-i 
Reggeons are increasingly important towards large z, and 
even dominate at z — > 1. 



B. Q 2 -dependence 

The full collections of ZEUS data named DIS with 
Q 2 > 2 GeV 2 can be binned in order to trace the Q 2 
dependence of the cross section. The result of such a 
binning is presented in Fig. [S] It demonstrate a clear in- 
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FIG. 9: (Color online) Q 2 dependence of the fractional DIS 
cross section of neutron production, integrated up to q™ ax = 
z x 0.69 GeV. ZEUS data [l| for Q 2 = 2.7, 8.9 and 40 GeV 2 
are shown by triangles, squares and round points respectively. 
The upper and bottom curves are calculated including all con- 
sidered mechanisms at Q 2 = 40 and 2.7 GeV 2 . 



crease of the fractional cross section with Q 2 , especially 
at z < 0.8. The variation with Q 2 of the calculated cross 
section originates solely from the absorption factor S 7 * 
(see Fig. 3]) which is increasingly important towards small 
z according to Eq. (|5Tj) . This naturally explains the ob- 
served trend of a weakened Q 2 dependence at large z. 



[601 ] . in the form of absolute values of the cross section. 
To compare with the ZEUS data and our calculations we 
normalized the HI data dividing by inclusive DIS cross 
section <7j nc = 112 pb [6l|. The results are plotted in 
Fig. [8] by squares. Data for the fractional cross section 
of both experiments agree with each other if integrate in 
the large interval qx < z x 0.69 GeV. However at small 
<7t < 0.2 GeV the HI cross section considerably exceeds 



C. q-r-dependence 

We also compare the qr dependence of the fractional 
differential cross section with few samples of ZEUS data 
presented in Fig. [101 Our results are shown by solid 
curves, which sum up the pion (long-dashed curves) and 
Reggeon (short-dashed curves) contributions. Appar- 
ently, the role of Reggeons increases with z, especially at 
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FIG. 10: Fractional differential cross section of leading neu- 
tron production as function of neutron transverse momen- 
tum at several values of z (solid curves). The contributions 
of pure pion exchange and of the iso- vector Reggeons, other 
than pion, are depicted by long and short-dashed curves re- 
spectively. ZEUS dara are plotted by round points. The large 
systematic errors (see [T3|) are not shown. 



large qx, where they significantly diminish the slope com- 
pared with the net pion contribution. Correspondingly, 
the calculated qx distribution acquires a bent shape, 
while the ZEUS data seem to prefer the Gaussian qx- 
shape. 

According to Eq. ([5]) at large z and small qx one ap- 
proaches the pion pole, l/(t + m%). Therefore, the qx- 
distribution of neutron production at large z is expected 
to be rather steep, since is controlled by the small pion 
mass. However, the relative contribution of natural par- 
ity Reggeons rises, and eventually, they take over at 
z — > 1. Therefore, one could expect a sudden drop of 
the slope of the (^-distribution. 

Few examples of the (^-dependence of the differen- 
tial cross section plotted by solid curves, are compared 
with ZEUS data in Fig. [TU] The contributions given by 
Eq. (f57|) and of Reggeons are plotted by long- and short- 
dashed curves respectively. Apparently, the Reggeons are 
important to achieve agreement with data. 

Data of ZEUS are also presented in 17[ as the slope 
of the qx distribution versus z, as is shown in Fig. 1111 
Theoretically, the slope is ill defined, if one does not 
specify in which interval of qx it was measured. The 
local slope, defined as 



da 



p—tn 



dzdq T 



(77) 



apparently may vary with qx- Although the results of 
the ZEUS experiment agree with q^-independent slope 
(within large systematic errors), this is certainly not the 




FIG. 11: (Color online) ZEUS data [17| for the <? T -slope of the 
differential DIS cross section of leading neutron production as 
function of z. The solid curves present theoretical predictions 
calculated with Eq. (73 at fixed q\ = 0.5, 0.1 and 0.2 GeV 2 . 



case in the theory. For this reason we do not perform an 
explicit comparison of our results with data, but present 
in Fig. [TT] the z-dependent slope B(z,qx) calculated at 
few fixed values of q\ = 0.5, 0.1 and 0.2 GeV 2 . Remark- 
ably, the results of calculations demonstrate flattening 
and a drop of B(z, qx) towards z = 1, similar to what is 
observed in data. Notice that these curves are not sup- 
posed to be directly compared with data, which cover 
an interval of qx dependent on z. The curves just show 
the uncertainty in the value of the slope, related to its ill 
definition. 



VI. DETERMINING THE PION STRUCTURE 
FUNCTION FROM DATA 

Now we are in a position to try to answer the question, 
whether the process of leading neutron production in DIS 
can be consider as a tool to measure the pion structure 
function at small x. The answer depends on the sensi- 
tivity of the cross section to the number of quarks in the 
proton at a soft scale, and on the involved theoretical 
uncertainties. 

If one trusted the way of calculations as is, then the 
number of quarks would affect only the ratio (|12p and 
the cross section as a simple rescaling factor. However, 
the procedure of interpolation of the absorption factor, 
Eq. (|51[) . between the known regimes of short and long 
coherence length, Eq. (|4"5)l . leaves some freedom in ad- 
justing the shape of the z-dependence of the fractional 
cross section. The parameter L is known only by the 
order of magnitude, it should be comparable with the 
nucleon size, i.e. L ~ 1 fm. So far we fixed it at this 
value, however, after rescaling the cross section assum- 
ing different number of quarks N%, we can readjust this 
parameter within a reasonable range in order to achieve 
a better agreement with data. 
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We considered two additional scenarios, which look to 
us extreme, = 3 and = 5, which correspond to no 
mesons, or in average to one meson in the light-cone wave 
function of the proton at a soft scale, respectively. Corre- 
spondingly, we adjusted the length scale at L = 2 fm and 
L = 0.5 fm respectively. The results a shown in Fig. [12] in 
comparison with data. Apparently, the solid curve does 
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• Even if the pion structure function is known, the 
fractional cross section of leading neutron produc- 
tion cannot be accurately predicted because of ini- 
tial/final state interactions generating absorptive 
corrections, whose magnitude is under debate. Al- 
though the projectile particle in 7*p — > Xn is a 
highly virtual photon, we expect the effects of ab- 
sorption at large z, to be nearly as strong, as in 
pp — > Xn. This happens due to formation of 
a strongly interacting colorOoctet dipole {qq\2 - 
{qq}%- The suppression factor S4 q (b) depicted in 
Fig. |3] has a magnitude similar to S5 q (b) in [2I [To|. 

• The coherence time Eq. (14T1) . which is the life- 
time of the strongly interacting projectile color- 
octet dipoles, is too short at medium values of z, 
and one should switch to the long-living qq fluctua- 
tions of the photon. The corresponding absorption 
factor S 7 * (b) is evaluated to be rather close to one, 
as is depicted in Fig. [4] The transition between the 
two regimes of absorption is illustrated in Fig. [5] 



FIG. 12: (Color online) Data are the same as in Fig. [8] upper 
panel. The dotted, solid and dashed curves are calculated 
with Ng — 3, 4, 5, and with L — 2, 1, and 0.5 fm, respec- 
tively. 



the best job describing the data, and being optimistic one 
may conclude that the version with N'P = 4 is preferable. 
However, being skeptical one may say that the results 
of calculations are too model dependent to make a solid 
conclusion. Moreover, the relation between the ratio of 
pion to proton structure function and the mean number 
of quarks at a soft scale is based on the model of two 
scales in the proton (24[, otherwise it may be broken. 



VII. SUMMARY 

To summarize, we highlight some of the results and 
observations. 

• Production of leading neutrons with fractional mo- 
mentum z — > 1 in DIS can serve as a way to mea- 
sure the pion structure function at small x. How- 
ever this method involves the several complications 
which are under investigation in this paper. 

• We expect a reduction of the ratio Fgfe)/!^ (x) 
compared with the usually assumed 2/3. The de- 
viation is due to the presence of a significant con- 
tribution of soft meson fluctuations in the proton 
wave function. While the contribution of the iso- 
vector mesons is constrained by the observed devi- 
ation from the Gottfried sum rule, the role of the 
iso-scalar mesons in the proton is less known. The 
calculations performed in this paper are done with 
the fixed ratio F£(x)/F%(x) = 1/2. 



• Other iso-vector meson exchanges, heavier than 
pion, are also important. The meson-nucleon cou- 
plings of natural parity p and ai Reggeons, which 
are predominantly spin-flip, were fixed by phe- 
nomenological Regge fits to high-energy hadronic 
data. The parameters of the unnatural parity a\ 
Reggeon, which is non-spin-flip, are not well con- 
strained by available data, and we fix them basing 
on the current algebra. Since the a\ pole contribu- 
tion was found to be very weak, we supplemented 
it by the p-it cut, and treated them together as an 
effective ai. Such an effective description was well 
tested in [2| with the data on neutron azimuthal 
asymmetry. The two sets of Reggeons have quite 
different intercepts and affect the neutron produc- 
tion cross section in different regions of z. Fig. [8] 
shows that p and 02 are important at z — > 1, while 
ai is large at smaller z. 

• Eventually, we additionally tested two extreme as- 
sumptions about the number of quarks in the pro- 
ton at a soft scale, Nfj = 3 and 5. In each case 
we readjusted the parameter L in Eq. (|5ip within a 
reasonable range, however, could not reach a good 
agreement with data, as is demonstrated in Fig.lT2"l 



Summarizing, our assumption Eq. (|21[) that the pion 
structure function at small x is twice smaller than the 
proton one, is well supported by the parameter-free cal- 
culations of absorptive corrections and contribution of 
Reggeons, providing a good description of ZEUS and HI 
data for leading neutron production in DIS, as function 
of z, Q 2 and qr- Nevertheless, we should admit that the 
test is not really precise because of many theoretical and 
experimental uncertainties involved into the calculations. 
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